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ABSTRACT 

We present results of 2D simulation of magnetorotational (MR) supernova accompa- 
nied by jet formation in the core collapse supernova explosion. Initial magnetic field 
used in the simulations has dipole-like symmetry. Contrary to the simulations of MR 
supernova with initial quadrupole-like magnetic field, where the matter was ejected 
mainly near the equatorial plane, in presence of the dipole-like initial magnetic field the 
supernova explosion is developing preferably along a rotational axis, and leads to for- 
mation of a protojet. We expect that protojet propagation through the envelope of the 
star will be accompanied by its coUimation. The magnetorotational instability (MRI) 
was found in simulations, similar to the earlier considered case of the quadrupole-like 
initial magnetic field. Our estimations show that the characteristic time for the re- 
connection of the magnetic field is much larger than the MRI development time. The 
supernova explosion energy for the dipole-like field is about 0.61 • lO^^erg, and about 
O.ISM© of mass was ejected during the explosion. 

Key words: MHD-supernovae: general. 



1 INTRODUCTION 



, The idea of the MR mechanism of the core collapse su- 
pernova, suggested bv iBisnovatvi-KogarJ lll97Cf) . was to get 
energy for the supernova explosion from the rotation of 
. the star due to magnetic field action. During the collapse 
' stellar rotation is becoming differential i.e. central parts 
of the star rotate much fa ster than the envelope. ID nu - 
merical simulat i ons m ade bv lBisnovatvi-Kogan et all il976l) , 
lArdelvan et al.l (|^7^ had shown that the differential rota- 
tion in presence of magnetic field leads to appearing of the 
toroidal component of the magnetic field, its linear growth, 
and following formation of MHD shock wave. In ID ap- 
proach the rotation was differential along r axis only (the 
star in ID was a cylinder) what means that the only Hr com- 
ponent of the magnetic field was winded up and produced 
the toroidal {H^p) component. In a 2D case the system has 
more degrees of freedom and the rotation is differential not 
only along the r axis. The differential rotation is using both 
components of the magnetic field {Hr, Hz) for the toroidal 
field amplification. 2D simulations of MR mechanism with 
the initial magnetic field of the quadrupol e -like type of sym- 
metry have shown ('Ardclia n et al.l (|200ff l , iMoiseenko et all 
(2004), Ardclian ct al. ( 200|)) that the amplification of the 
toroidal component of the magnetic field leads to formation 
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of a compression wave, which moves along steeply decreas- 
ing density profile and quickly transforms into the fast MHD 
shock wave. Due to the quadrupole-like initial magnetic field 
the supernova explosion was developed predominantly along 
the equatorial plane. The explosion energy for that case was 
about 0.6- 10^^ erg. The value of the ejected energy is enough 
for the explanation of the core collapse supe rnova explosion. 
The m ain qualit ative difference betwee n ID jArdelvan et alJ 
(|l973)) and 2D iArdelian et al.l i200^ ) MR supernova sim- 
ulation is an appearance of the magnetorotational instability 
(MRI) which significantly reduces the time of the growth of 
the toroidal magnetic field up to MR supernova explosion. In 
ID case the MRI is not developed due to the small number 
of degrees of freedom. 

The first simulations of the MR p rocesses in stars have 
been done bv lLeBlanc fc Wilson|[l97ol after which MR pro- 
cesses in the stars in relation to the core collapse supernova 
explosion had been simulated by 'Bisnovatvi-KoEane^^ 
d976, Ardclyan c t al. 1979, Miiller & Hillcbrandt Il97gl . 
lohnishilligsA ISvmbalistvlll984l Recently the interest to the 
MR processes (especially in application to the core collapse 
supernova) was recommenced due to increasing number of 
observational data about asymmetry of the explosion, and 
possible colli mated ejecta in connec ti on with cosmic gamma 
ray bursts jArdeban et alJ l200d: lAkivama et alj l2003l: 
l Yamada fc Sawaill2004 iTakiwaki et alJl2004l: iKotake et alJ 
I2OO4I: lArdeban et alJl2005l: lYamasaki fc Yamadd l2005^ . In 
the paper bv lBignova^IJ<ogai^^Ti^uko^T2OT it was es- 
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timated that a significant part of the type Ib,c supernovae 
could explode as magnetorotational supernovae. 

In the present paper we describe the results of the 2D 
simulations of the MR supernova explosion with the ini- 
tial dipole-like magnetic field. The qualitative difference be- 
tween our simulation of the same problem with the initial 
quadrupole-like magnetic field is the formation of mildly col- 
limated prot ojet, directed alo ng the rotational axis. As in 
the paper bv lArdelian et all l2005i) we have found a devel- 
opment of the MRI leading to rapid explosion even at a rela- 
tively weak initial magnetic field. In the simulations we have 
used an ideal MHD with the infinite conductivity. In spite 
of the fact that we supposed the conductivity to be infinite, 
in reality it is very big, but finite. That means that during 
winding up the force lines of the magnetic field and the MRI 
development, the reconnection of the magnetic field lines is 
possible. If the reconnection time would be smaller or com- 
parable to the characteristic time of the MRI development 
then the reconnection process could reduce effectiveness of 
the MR supernova mechanism. 

We have done estimations for the reconnection time and 
have found that it is much larger than the time of the de- 
velopment of the MRI, and will not influence significantly 
on the MR supernova explosion (MRE). MRI in app l icatio n 
to the core collapse su pernova was studied bv lSoruid J2002D . 
lAkivama et al.l jioO^)- 

For the simulations we have used implicit Lagrangian 
numerical method on triangular grid of variable structure. 
This method was used for the simulations of a number of 
different astrophysical problems: simulations of the collapse 
of co ld rapidly rotating protostellar clouds ( Ardcli an et alJ 
il996) '). magnetorotational processes for the colla psing mag- 
netized protostellar cloud (lArdelian et al.l ll200(J)'l. core col- 
lapse and formation of the neutron star i^rdelian et alJ 
l|200^), MR super nova simulations with q uadrupole-like ini- 
tial magnetic field llArdelian et alJ tOQ^) ]. 



2 FORMULATION OF THE PROBLEM 

2.1 Equations of state and neutrino losses 

For the simulations w e use the equation of state from 
lArdelvan et"al] digSTaTl : 



P = P{p,T)= Po(p) + pUT+^, 



3(1) 
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Here 5R is the gas constant taken equal to 0.83 x 
10* cm'^s^'^A', a is the constant of the radiation density, 
P is pressure, p is density, and T is temperature. In the 



expression -Po(p) the value p was identified with the total 
mass-energy density. For the cold degenerated matter the 
expression for Po( p) is the ap proximation of the tables from 
iBavm et al.l lll971l) : iMalone et al. (197^. 

In the neighborhood of points p = pfc in equation Q 
the function P o (p) wa s smoothed in the same way as in 
lArdelvan et"ai] 1(19873), to make continuous the derivative 
dPo/dp : 



Mp) 



where 
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The specific energy (per mass unit) was defined ther- 
modynamically as: 



e = eo(p) + I^T + — + eF.(p,T). 
2 p 

The value eo(p) is defined by the relation 
' ^'o(p). 



eo(p) 



idp. 



(3) 



(4) 



The term epe in equation @ is responsible for iron 
dissociation. It is used in the following form: 



£Fe(p,T) 



T-Tn 



A nip \ TiFe — Toj 



(5) 



It is supposed that in the region of the iron dissocia- 
tion the iron amount is about 50 per cent of the mass, 
Eb,Fe = 8 X 10~"'erg is the iron binding energy, ^ = 56 is the 
iron atomic weight, and rUp — 1.67 x 10~^*g is the proton 
mass, ToFe = 0.9 x 10^°K, TiPe = 1.1 x 10^°K. For the nu- 
merical calculations formula @ has been slightly modified 
(smoothed): 



eFe{p,T) 



Eh 



A m- 



p 2 L 



T — ToFe 
TlFe — Toi 



}-(6) 



The neutrino losses for Urea processe s are u sed in the 
form, taken from iBisnovatvi-Kog an et alJ (|l976l ). approxi- 
mating the table of l^emovE^^^ ^1969^ : 



/(p,r) 



1.3 • io^ae(r)r 

+ (7.1 • 10-5pT^)§ 



erg - g 



ae(T) = { 664.31 



(7) 



(8) 



1, T< 7, 

51.024(r - 20), 7<T<20, 

664.31, T > 20, 

T = T- 10"^ 



The neutrino losses from pair annihilation, photo pro- 
duction, and plasma were also taken into account. These 
types of the neutrino losses have been approximat ed by the 
interpolation formulae from ISchindler et alJ lll987f) : 



Qtot — Qpair Qphoto ~\~ Qplas 



(9) 
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The three terms in equation @ can be written in the fol- 
lowing general form: 



fflo + aiC + "2^ 



(10) 



+ bia + b2a^ + baa^ 

For d — pair, K{p, a) = g{a)e~'^°', 

, . , 13.04 133.5 1534 918.6 
g(a) = \ — + — — H g- + — g-; 

Q,4 Q,b Q,» 

For d = photo, K{p, a) — {p/ fiz)a~^; 
For d = plasm, K{p, a) — (p/pz)^; 

Here, pz = 2 is the number of nucleons per electron. 
Coefficients c, Uj, and bj for th e different losses are given in 
the Table m of FSchindler et"al] (|T987). The general formula 
for the neutrino losses in a nontransparont star has been 
written in the form, used bv lXrdelian et al. (20041 : 



F{p,T) = {f{p,T) + Qtot)e-To. 



(11) 



The multiplier e"!^ in equation (IllL where = S^nly, 
restricts the neutrino flux for non zero depth to neutrino 
interaction with matter Tv The cross-section for this inter- 
action Sv was presented in the form: 

" (0.5965 ■ 1010)2 ' 
and the concentration of nucleons is 

nip 

The characteristic length-scale Iv, which defines the depth 
for the neutrino absorbtion, was taken to be equal to the 
characteristic length of the density variation as: 

P ^ P_ 

Vp\ 



I. 



(12) 



{{dp/drY + (dpjdzff^'^' 

The value 1^ monotonically decreases when moving to the 
outward boundary; its maximum is in the centre. It approxi- 
mately determines the depth of the neutrino-absorbing mat- 
ter. The multiplier 1/10 in the expression e"^"!^ was ap- 
plied because in the degenerate matter of the hot neutron 
star only some of the nucleons with the energy near Fermi 
boundary, approximately 1/10, take part in the neutrino 
processes. 

2.2 Basic equations 

Consider a set of magnetohydrodynamical equations with 
self-gravitation and infinite conductivity: 

dx 



At 



V, 



dp 

At 



' pV • V = 0, 



H-H\ V-(H®H) 



dv ,/„ H-H\ V- 



47r 



pV^, 



(7) = H • Vv, A$ = 47rGp, (13) 

P^ + ^'V-v + pF(p,T) =0, 
P^P{p,T),e^e{p,T). 



here %i +v- V is the total time derivative, x = (r, i^j, z). 



V — {vr,v^,V2) is the velocity vector, p is the density, P is 
the pressure, H = {Hr, H<f,, Hz) is the magnetic field vector, 
<E> is the gravitational potential, e is the internal energy, G 
is gravitational constant, H(8)H is the tensor of rank 2, and 
F{p, T) is the rate of neutrino losses. 

r, ip, and z are spatial Lagrangian coordinates, i.e. 
r = r(ro,((3o, and ZQ,t), ip = ip{ro, ipo, zo,t), and z = 
z{ro, ifiOjZct), where ro,(fio, zo are the initial coordinates of 
material points of the matter. 

Taking into account symmetry assumptions (^ =0), 
the divergency of the tensor H (g) H can be presented in the 
following form: 



1 a(rHr-Hr-) I a(H^Hr) 
Or dz 



V.(H®H)= l2iLMLlM + 2iE^ + lH^H.^ 



a. 

1 ajrHrH^) 



ar 



+ 



dz 



Axial symmetry 



— 0) and symmetry to the equa- 
torial plane are assumed. The problem is solved in the re- 
stricted domain. At f = the domain is restricted by the 
rotational axis r ^ 0, equatorial plane z ^ 0, and the outer 
boundary of the star where the density of the matter is zero, 
while poloidal components of the magnetic field Hr , and Hz 
can be non-zero. 

At the rotational axis (r = 0) the following boundary 
conditions are defined: (V4')r = 0, Vr = 0. At the equato- 
rial plane {z — 0) the boundary conditions are: (V'I>)z = 
0, Vz = 0. At the outer boundary (boundary with vacuum) 
the following condition is defined: Pq^^^j- boundary ~ 

We avoid explicit calculations of the function eo(p) in 
equation (01, because this term is eliminated from (I13II due 
to adiabatic equality: 



deo fo dp 

p— = — =PoV- V, 

at p at 



(14) 



determining the fully degenerate part of the equation of 
state. Therefore, we define 

3 ctT'^ 

e* = -3fJr+ +eFe{p,T), 

2 p 



P* = pKT + 



The equation for the internal energy in equation (1131 can be 
written in the following form: 

P^+P*V-v + pF(p,r) =0. (15) 



2.3 Initial magnetic field 

The initial poloidal magnetic field is defined sim ilar to 
our pr evious papers (Ardelian et al. (2005), Ardclia n et alJ 
(|2000|)) by the toroidal current j'^ using Bio-Savara law. The 
toroidal current which determines the initial magnetic field 
should be defined in the upper and in the lower hemispheres. 
The field with the quadrupole-like symmetry is formed by 
the toroidal current antisymmetrical to the equatorial plane. 
The dipole-like magnetic field is formed by current symmet- 
rical to the equatorial plane. To define the initial poloidal 
magnetic field we have used the following toroidal current 
(in nondimensional variables: 
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Table 1. 


Coefficients for equation 


1101 fromlSchindler et al 


.Nl98'7t). 
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ai 
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10** X ^ T ^ 10^" K 


pair 
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5.026(19) 
3.897(10) 
2.146(-7) 
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1.568(21) 1.2383 
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1.653(-8) 2.581(-2) 


-8.1141(-1) 

7.483(-3) 
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0.0 

3.061(-4) 
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0.56457 




Figure 1. The initial poloidal dipole-like magnetic field configu- 
ration. 



where 



ji^forz >,< 0; 



for 



(r -0.15)"+ 0.2(^)2 <:0.3^ 



Aj X cos 



(r--0.15)" + 0.2(z)^ >0.3" 



(r - 0.15)2 + 0.2«2 



0.32 



■ d -u 

Jip J:p ■ 



The coefficient Aj serves for adjusting the initial magnetic 
energy. In the dipole-like type of symmetry 1161 there is 
Hr = at the equatorial plane {z — 0) (Fig. • It is im- 
portant to have a force-free or a balanced initial magnetic 
field. Using of the nonbalanced initial magnetic field leads 
to appearance of artificial flows due to "turning on" effect. 
We start the simulation of the MagnetoRotational Explo- 
sion (MRE) from the "equilibrium" state to avoid numeri- 
cal problems, arising from the nonstationarity, especially in 
outer region of the core. These regions give a small input 
into the total energy balance of the explosion. 

To construct a configuration, where magnetic forces are 
balanced with other forces (i.e. gravitational force, gradient 
of ga s pressure and centrifu gal fo rce) we follow the proc edure 
from lArdelian et al.l tOOdti and lArdelian et alJ (l2005l) . The 



Figure 2. Time evolution of the toroidal (solid line) and poloidal 
(dashed line) parts of magnetic energy during MR explosion with 
dipole-like initial magnetic field. 



ratio between the initial magnetic an d gravitational energie s 
was chosen to be equal to 10~^ , like in lArdelian et all i2005l) , 
for better comparison of the present results with the results 
of simulations of MR explosion with the initial quadrupole- 
like magnetic field. The initial poloidal magnetic field in the 
center of the star at start of the evolution of the toroidal 
magnetic field is ~ 3.2 x 10^^ G. 



2.4 Core collapse simulation 

We first have done simulation of the rotating nonmagnetized 
core collapse llArdelian et al.ll2004ft . Initially the ratios be- 
tween the rotational and gravitational energies and between 
the internal and gravitational energies of the star had been 
chosen as: 



Er 



Ea 



0.0057, 



Ei, 



Ea 



0.727. 



The rotating core collapse resulted in formation of the 
proto neutron star with an extended envelope. The bounce 
shock which was formed during the collaps e, goes through 
the envelope. At the end of calculations in lArdelian et alJ 
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l|2004|) at t = 0.2565 s the bounce shock reaches the outer 
boundary of our computational domain. The shock leads to 
the ejection of 0.041 per cent of the core mass and 0.0012 per 
cent (2.960 x 10** erg) of the gravitational energy of the star. 
The amounts of the ejected mass and energy are too small 
to explain the core collapse supernova explosion. After the 
core collapse {t — 0.261 s) we obtain a differentially rotating 
configuration. 

As in previous simulation of the MRE we divide the 
whole process into three separate parts. The first part is 
the simulation of the core collapse. The second part is the 
construction of the balanced star with the dipole-like mag- 
netic field. The third part is the calculation of MRE of this 
configuration. 

The initi al model for simulation of MRE was taken from 
the paper of lArdelian et all ()2004h . where the first part of 
the process (collapse) was calculated. It is represented by 
differentially rotating hot proto neutron with central den- 
sity ~ 2 X lO^'* g/cm^, central temperature 3.5 x 10^" K, 
and angular velocity near the center, corresponding to the 
period P ~ 10~^s. At equatorial radius ~ 10 km, contain- 
ing ~ 18% of the mass the rotational period was equal to 
1.5 X 10~^ s, and increased up to ~ 35 s at the outer bound- 
ary at the equator. T he collapse calculations in the paper of 
lArdelian et alJ (|2004 ) have been started from rigidly rotat- 
ing iron "white dwarf star with the mass ~ I.2M0, 20% 
larger than Chandrasekhar limit for the pure iron star, with 
the rotational period P ~ 0.4 s. The outer boundary is ex- 
panding during the collapse, leading to large increase of the 
rotational period. 

After the core collapse and formation of a steady state 
differentially rotating configuration, the initial poloidal mag- 
netic field, defined by th e toroidal curre nt 116( was 'turned 
on'. As in the paper bv lArdelian et al at the sec- 

ond step the toroidal component was switched off, to obtain 
a balanced magnetized configuration, and the evolution of 
only poloidal components {Hr, Hz) of the magnetic field was 
permitted. 



3 NUMERICAL METHOD 

The numerical method used in present simulations is based 
on the implicit operator-difference completely conservative 
scheme on the Lagrangian triangular grid of variable struc- 
ture. The implicitness of the applied numerical scheme al- 
lows to make large time steps. It is important to use implicit 
scheme in such kind of problems due to the presence of two 
strongly different timescales. The small timescale is defined 
by the huge sound velocity in the central parts of the star. 
The big time scale is defined by the characteristic timescale 
of the evolution of the magnetic field. Conservation prop- 
erties of the numerical scheme are important for the exact 
fulfillment of the energy balance and divergence-free prop- 
erty of the magnetic field. 

During the simulation of the MRE the time step for the 
implicit scheme was ~ 10 300 times larger than time step 
for the explicit scheme (CFL stability condition). It means 
that the total number of time steps is 10-^300 times less than 
it would have been done for explicit scheme, what allows 
us to decrease time approximation error. We did not make 
direct comparison of CPU time per time step for our im- 



plicit scheme and an explicit scheme. Our estimations show 
that the total number of arithmetic operations for the im- 
plicit scheme is ~ 20 times larger than it is required for 
explicit scheme. Considerable decrease of the required num- 
ber of time steps leads to corresponding reduce of the time 
approximation error. 

Grid reconstruction procedure applied here for the re- 
construction of the triangular lagrangian grid is used both 
for the correction of the "quality" of the grid and for the 
dynamical adaptation of the grid. 

The method applied here was developed, and 
its stability was inves t igated in the pape r s by 
lArdelian fc Kosmachevskiil lll995l) . lArdelian et alJ lll987t) 
and reference s therein. It was te sted thoroughly with 
different tests (|Ardelian et alJ (|200(t '). 

4 RESULTS 

4.1 Magnetorotational supernova explosion and 
protojet formation 

After the core collapse we get a steady differentially rotating 
configuration which properties are described in the section 
12.41 The ratio of the rotational energy and gravitational en- 
ergy is equal to 0.073. The star remains in such a state pretty 
long time. We have done about 10000 time steps (it corre- 
sponds to 0.03 s of physical time) and did not notice any 
significant changes of the parameters of the star. While the 
inclusion of even weak initial poloidal magnetic field leads 
to drastic changes. 

After formation of the balanced configuration with the 
poloidal dipole-like magnetic field, at the moment of 'switch- 
ing on' the equation for the evolution of the toroidal mag- 
netic field we start counting the time anew. 

The toroidal magnetic field component appearing due to 
the differential rotation is amplifying with time. The fastest 
amplification of the Hip takes place in the regions of the max- 
imum of the value rH -grad^Vip/r) in the equation of the evo- 
lution of Hip. At the initial stage Hp grows linearly with time 
(correspondingly toroidal magnetic energy Hp / (Svr) grows 
as a quadratic function) (Figl^J. Due to appearance of the 
MRI (which will be described in the following section), af- 
ter the stage of a linear growth, the toroidal magnetic field 
starts to grow much faster (exponentially) . At this stage the 
poloidal magnetic field also begins the exponential growth. 

The toroidal magnetic energy reaches its maximal value 
^ 3.2 X 10^" erg at t = 0.07 s. This value is « 66% of 
the maximal value of the toroidal energy in the case of ini- 
tial quadrupole-like m agnetic field, where th is maximum was 
reached at t = 0.12 s (lArdelian et alJ l|200^). The maximal 
values of the toroidal mag netic field Hip « 1.8-10^^ G (in the 
case of the quadrupole this maximal value was ~ 2.5 x 10^^ 
G) are reached at the distance ~ 10 km from the center of 
the star. 

After reaching its maximal value the toroidal magnetic 
energy decreases with the time. At the developed stage the 
poloidal magnetic energy reaches ~ 4.5 x 10^° erg, and does 
not show any substantial decrease with time. 

At t = 0.07s a compression wave appears moving from 
the central parts of the star due to increasing of the magnetic 
pressure. Moving along a steep density profile the compres- 
sion wave transforms into the fast MHD shock wave. The 
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Figure 3. Time evolution of the rotational and radial kinetic 
energy during MR explosion. 



Figure 4. Zoomed part of the time volution of the radial kinetic 
energy (Fig. |3J during MR explosion. 



development of the MHD shock outwards leads to the ejec- 
tion of mass and energy. 

During the MR explosion a part of the rotational en- 
ergy (« 10%) is transforming into the radial kinetic energy 
(which could be related to supernova explosion energy). The 
neutrino transport effects are taken into account approx- 
imately in flux limited approximation variant (see section 
12. H . A larger fraction of the rotational energy is lost by 
neutrino emission, while most part of the angular momen- 
tum is carried away by the ejected matter. 

At the Fig|3| the time evolution of the rotational and 
radial (poloidal) energies is represented. At the final stage 
of the MRE simulations, the excitation of eigenmodes of the 
proto neutron star happens. The plot of the radial kinetic 
energy (Fig. 2J shows that this eigenmode has a period of 
~ 1 ms. 

At the Fig. |3 the kinetic energy of the core is repre- 
sented. This plot was made by subtracting the kinetic en- 
ergy of the matter expanding by the supernova shock from 
the total kinetic energy of the star. The amplitude of the 
kinetic energy of the core osci l lation s is ~ 1.3 • 10*** ergs. 
In the paper bv lBurrows et aP i2006l) a mechanism of core- 
collapse supernova explosions was suggested, based on the 
acoustic energy extracted from the core. 

Our simulations have shown that the amount of the 
acoustic energy generated in the collapsing core is much less 
than the expected energy of core collapse supernova, note 
that our mathematical method was different, and we have 
used more simplified desc ription of the physical processes 
than lBurrows et alj J2006t) 

At the Figs|Sl |7|the time evolution of the neutrino lu- 
minosity and neutrino losses during MR explosion are rep- 
resented. 

The amounts of ejected mass (r; O.ISMq) and energy 
(6.1 X lO^^erg) are close to the corresponding values in the 
case of MR supernova explosion with the initial quadrupole- 
like magnetic fleld. The plots of the time evolution of the 
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Figure 5. The time volution of the kinetic energy of the core 
only during MR explosion. 



ejected mass and energy for the dipole field are represented 
at FiglHland Figl^l 

The time evolution of the velocity field and temperature 
are represented in Figs llUl and llll At the upper plot of the 
Fig llUI a zoom of the central part of the velocity field is 
represented. The shape of the supernova MHD shock is not 
spherical. Near z axis a protojet is forming. At the developed 
stages of MR explosion the shape of the shock front becomes 
more spherical, while the plot for the time evolution of the 
specific angular momentum Fig ll2l shows that the angular 
momentum is extracted mainly along the axis of rotation z. 

A development of the supernova shock leads to forma- 




time, s 



Figure 6. Time dependence of the neutrino luminosity 

rA'I 

core p(^p 'j''j^jyi during MR explosion. 
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Figure 8. Time evolution of the ejected mass during MR explo- 
sion. 




time, s 



Figure 7. Time dependence of the neutrino losses during MR 
explosion. 




time, s 



Figure 9. Time evolution of the ejected energy during MR ex- 
plosion. 

more complicated what means that the resulting shapes of 
the MR supernova explosion could be rather different. 



tion of a mildly collimated protojet. This protojet could be 
coUimated stronger during passing through the vast enve- 
lope of the star and formation of the funnel flow. 

As we have shown in the paper lArdelian et alJ j^QQ^ 
the MR mechanism with the initial quadrupole-like mag- 
netic field leads to the explosion which develops prefer- 
ably near the equatorial plane. The shape of the MR ex- 
plosion qualitatively depends on the initial configuration of 
the magnetic field. The magnetic field with a pure dipole or 
quadrupole types of symmetry in the star is a simplification. 
In reality the structure of the magnetic field could be much 



4.2 Appearance of the magnetorotational 
instability (MRI) in 2-D picture 

It was shown in ID si mulations of MR supernova explosion 
lArdelvan et al.l (^7^ that the generated toroidal magnetic 
field is growing linearly with the time. Our recent 2D MR 
supernova s imulations with the qu adrupole-like initial mag- 
netic field jArdelian et al.l i2005ft l have shown that, after 
about 100 rotational periods of the central core, its linear 
growth changes to the exponential one, and the poloidal 
components also start to grow exponentially. The reason for 
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Figure 10. Time evolution of the velocity field for time moments 
t = 0.075s, 0.1s, 0.25s 
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Figure 12. Time evolution of the specific angular momentum 
: 0.075s, 0.1s, 0.25s. 



Figure 13. The MRI dcvclopinciit for time moments t = 
0.0045s, 0.018s, 0.042s. Gray scale is the toroidal field Hep levels. 
Arrows show a direction and strength of the poloidal magnetic 
field Hr, H^. 
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that exponential growth i s an onset of the variant of MRI, 
which was investigated bv lTavleJ (^7^. The qualitative ex- 
planation of MRI development for the MR supernova in 2D 
case was given by Ardelian ct al. ( 2005). Appearance of the 
MRI significantly reduces the evolution time of the MR su- 
pernova explosion. The picture of the MRI development in 
the present simulations is represented in Fig |13l 

Let us define a as a ratio between the magnetic EmagO 
and gravitational Egravo energies at the moment of "turn- 
ing on" of the magnetic field. It was found in ID simula- 
tions that the time of the field amplification until the MR 
explosion, which is equal approximately to the time of the 
whole explosion t^^pi, depends on a as texpi ~ The 
results of 2-D MR supernova simulations with the initial 
quadrup ole-like magnetic field ha ve been represented in the 
paper bv lMoiseenko et alj (1200411 for a wide range of varia- 
tion of a = 10"^ 10"^^. It was found, that for a <~ 10"" 
the growth of texpi with decreasing a is becoming much 
slower, and may be well approximated by a logarithmic for- 
mula t^xpi ~ I log a I for small q ^ 1. In the case of the 
dipole-like initial magnetic field we have made simulations 
for the same range of variation of the a parameter, and have 
found a similar dependence which is represented in Fig |14l 

The logarithmic functional dependence at very small a 
has a simple qualitative explanation. First, we should have 
in mind, that the development of MRI starts at t = tmri, 
approximately at the same moment for all initial values of 
the magnetic field. The MRI development begins in the part 
of the star, where the ratio of the toroidal and poloidal mag- 
netic fields reaches a definite value / ~ few tens, what hap- 
pens after about 100 rotational periods of the inner core. In 
another words, the time of the field growth until the begin- 
ning of MRI tmri does not depend on a, if a is small enough. 
At larger a the explosion starts before the beginning of MRI. 

Second, we should take into account that MR explosion 
happens approximately at the same ratio F < 1 of the mag- 
netic Emag and internal Eint (gravitational Egrav) energies 
of the star, also for all a, EmagE ~ FEgrav The magnetic 
energy Emag is growing exponentially with time during MRI 
development at t > tmri, 



E„ 



= exp[7,„(t — tmri)], 



(17) 



fE 

where 7^ is an increment of the MRI instability. Therefore, 
the time between the the start of MRI, until the moment of 
the explosion texpi , is growing logarithmically with decreas- 
ing a 



, , 1 1 EuiagE 

f^expl <^mri — iOg 

J -tymagO 



log-TTT^ 

'^m J i-^magO 



1 F EqravO 

— log 



1 n ^ 1 ^ 
— (log -J - log a). 

7 m / 



(18) 



Here it was taken into account that a relative change of the 
gravitational energy of the star during MagnetoRotational 
Explosion (MRE) is rather small, and Egrav{t) ~ Egravo- 
At larger a the MRE happens before MRI instability be- 
gins, texpi < tmri, with t^^pi ~ l/y'Tx, and for small a, when 
tsxpi >> tmri, the exploslou time is growing logarithmi- 
cally. This qualitative explanation is valid if the explosion 
happens during a linear stage of the MRI development. On 
this stage the magnetic field is growing exponentially with 




Figure 14. The dependence of the time evolution of the MR 
explosion on a = f?"'"" for the case of the initial dipole-like 
magnetic field. 



the increment roughly defined by the values of the config- 
uration at the beginning of instability. When the magnetic 
field energy roughly approaches the rotational energy of the 
star, its growth is stopped, resulting in nonlinear saturation 
of MRI. Our numerical simulations have shown, that the 
explosion occurs before the nonlinear saturation of MRI. 



4.3 Influence of numerical dissipation, and 
magnetic field reconnection 

We have done MR supernova explosion simulations with the 
initial dipole-like magnetic field (Fig0 at a different num- 
ber of knots in our triangular grid, starting from 1500 up to 
18000 knots. We have found that the moment of the onset 
of the MRI depends visibly on the number of grid points. 
The less knots we have in the grid, the longer time it takes 
for starting the MRI. We have found that the results of the 
simulations (namely the time of the development of the MR 
supernova explosion) do not change when the total number 
of the knots Nn exceeds 15000. The results of the simula- 
tions described here are obtained with the grid containing 
about 15000 knots. The dependence (converging) of the time 
of the explosion t^^pi on the number of knots has a pure 
numerical origin, and is connected with a numerical dissi- 
pation. Such dissipation has a stabilizing influence on the 
MRI, and its onset is shifted to larger values of / - the ratio 
of the toroidal and poloidal magnetic fields. Therefore, the 
total time of MRE is increasing with decreasing of the num- 
ber of knots, due to increasing of the time tmri- Numerical 
dissipation influence also on the field growth during MRI 
development, decreasing the increment 7™,, giving the same 
effect. As follows from our simulations the influence of the 
numerical viscosity becomes unimportant at N„ ^ 15000, 
when the time of viscous (numerical) dissipation starts to 
exceed all characteristic times of MR explosion. 

In the simulations we supposed that the matter of the 
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star has infinite conductivity, although in reahty the conduc- 
tivity is finite. In the case of a finite conductivity and chaotic 
structure of the magnetic field the reconnection of the mag- 
netic field could be important, leading to a real physical 
dissipation. 

The theory of the magnetic field reconnection is not 
complete up to now. To analyze the importance of the mag- 
netic field reconnection in MRE, let us estimate (at least 
roughly) a characteristic time for this process. The growth 
of the magnetic field is accompanied by appearance of the 
chaotic small-scale structure of the magnetic field. If the 
characteristic reconnection time in smaller scales is compa- 
rable with characteristic time ~ l/7m of the magnetic field 
amplification due to MRI action, then reconnection could 
lead to a delay in the magnetic field growth, and hence could 
decrease the efficiency of MR explosion. 

For the estimation of the characteristic time for the 
magnetic field reconnection we use the formula for the re- 
connection time of the quickest - Petsch e k-typ e reconnection 
JPriest fc SowardI Jl97g) . lGaleev et al.l t97^ ): 

_ 4(ln(Re„0 + 0.74) 

'reconn — t i 5 x^"^} 

tvvaI~ 

where va = H /(A-KpY^'^ is the Alfven speed, Rcm = 
A-K(jul/(? is the magnetic Reynolds number, I is the length. 
It is the velocity, p is the density, a is the microscopic 
plasma conductivity, c is the light speed. In the region 
where the MRI is developing we have the following param- 
eters: M = 0.2 • 10^° cm/s, I = lO^cm, p = 7 • lO^^g/cm^, 
H = 0.17- 10^*0. 

The plasma c onductivity due to binary collisions a is 
defined as (see e.g. lKrall fc Trivelpiecd il973t) ') 



/2aeT\3/2 



(16V¥)Ze21nA V J 



(20) 



where rUe = 9.1 ■ 10~ g is the electron mass, ee— 1.38 ■ 
10~^^erg K"-^ is the Boltzmann constant, e = 4.8 ■ 10"^" is 
the electron charge, Z = 27 is the ion charge, InA — 10 is 
the Coulomb logarithm, T — 10^^ K is the temperature. The 
above expression gives the conductivity of electrons in a non- 
relativistic, non-degenerate plasma. Relativity may be taken 
approximately by substituting the light speed instead of 
the expression in brackets in I20L and increasing the elec- 
tron mass in the numerator by relativistic factor 7 ~ 20. 
Substituting all these data in formula 1201 we find that con- 
ductivity is (7 ~ 8 ■ lO^^s"^. Note, the degeneracy of the 
electrons increase its conductivity, so this estimation may 
be considered as a lower limit for a. The magnetic Reynolds 
number in our case is Rcm ~ 10^^. Substituting these data 
into the formula 11911 we find the absolute lower limit for the 
characteristic reconnection time as 



' reconn 



5 s 



Our numerical simulations show that the characteristic time 
of the development of MRI is much smaller (~ 10 times) 
than the characteristic reconnection time. We can conclude 
that the reconnection processes will not suppress magnetic 
field amplification and hence MR supernova explosion. The 
formula 1191 for the characteristic time of the reconnection 
used here for the Petschek-type reconnection is the lower 
limit for the reconnection time. We can not assert that this 



type of the reconnection will be realized in the case of MR 
supernova mechanism while in the case of the realization of 
any other reconnection model the reconnection time will be 
even larger. 



5 DISCUSSION 

The results of the 2D numerical simulations of the MR su- 
pernova explosion with the initial dipole-like magnetic field 
have shown that MR supernova explosion is sensitive to the 
magnetic field configuration. MR supernova explosion with 
the initial quadrupole-like magnetic field develops mainly 
near the equatorial plane. In the case of the dipole-like mag- 
netic significant part of the ejected matter obtains a velocity 
along the rotational axis. The total energy of the explosion 
in the MR mechanism does not depend significantly on the 
topology of the initial magnetic field. The explosion energy 
for quadrupol e-like magnetic field is approximately equal to 
0.61 • 10^^ erg lArdelian et alJ ^QQ^, and for the dipole-like 
magnetic field its value is about ^ejected ~ 0.5 • lO^'^ erg. 
The amount of the ejected mass is approximately the same in 
both cases Mejected ~ O.14M0. Comparison of the explosion 
times for the dipole-like and the quadrupole-like fields of the 
same initial magnetic energy, shows that in the quadrupole 
case the explosion is developing faster. The time of the ex- 
plosion for the a — 10~® is about ~ 0.12 s in the dipole case, 
and is about ~ 0.06 s in the case of the quadrupole. 

In both cases we have done 2D simulations in a wide 
range of the initial magnetic field strength. The parameter 
a (the ratio between the initial magnetic energy and the 
gravitational energy of the star at the moment of "turning 
on" of the magnetic field) was chosen as 



EniagO 



10" 



10" 



(21) 



Comparison o f the results of ID simulations from 
lArdelvan et al.l lll979ll . with a dependence of the MR ex- 
plosion time on a as 

with corresponding 2D results ("Fig |141 . shows a qualitative 
difference between them. The reason is a development of the 
MRI in 2D-case, which reduces drastically the MR explosion 
time texpi- For the values of a >~ 10"'' the dependence 
1221 holds approximately, while for smaller a values the MR 
explosion time depends on a in the following way: 



expl 



1 log a\ 



(23) 



The axial protojet forming in our simulations is not narrow, 
while propagating through the envelope of the massive star 
this protojet could be collimated, or suffer mass entrainment 
and deceleration. 

The collapse does not lead to strong amplification of the 
toroidal field, because its time in the central parts is com- 
parable with rotation time. Besides, the initial configuration 
is probably stable against MHD instabilities and (in reality) 
contains both parts, toroidal and poloidal, of the magnetic 
field. The relation between these components, and details 
of the initial angular velocity distribution (uniform in our 
case) may lead to different situations, when MRI starts to 
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develop during the collapse, or after many rotational periods 
of the new born neutron star . Here we investigate the second 
case, as the initial step to the problem. The time separation 
of the hydrodynamic collapse from MRI and MRE helped 
us to understand the physical picture of MRI development, 
and to obtain the estimations for the energy production in 
MRE. 

We plan to extend our calculations to a more realistic 
case when magnetic field will be included at the very begin- 
ning of the collapse. We expect to obtain more complicated 
picture of magnetohydrodynamical processes, which should 
not infiuence strongly on the energy production, depending 
mainly on the initial rotational energy and angular velocity 
distribution. 

In the paper bv lSawai et al.l ll2005ll a similar problem 
with rather strong initial magnetic field was simulated, and, 
contrary to our results, the authors did not find the de- 
velopment of MRI. Their simulations of the supernova ex- 
plosion with weaker magnetic field did not lead to the de- 
velopment of MRI either and hence they did not get super- 
nova explosion for a weak initial magnetic field. The authors 
iSaw ai et al. 2005) claim that MRI did not develop in their 
simulations due to the insufficient spatial resolution. The 
absence of the MRI in their calculations most probably is 
connected with a rather large numerical viscosity of the nu- 
merical scheme they used. The numerical viscosity of the 
simulations depends on characteristic size of the mesh. The 
bigger is the size of the cell the larger is the numerical vis- 
cosity. In the section ll^ we have discussed the change of the 
MRI appearing time moment in dependence of the spatial 
resolution of our grid. In the case when the spatial resolu- 
tion of the grid is very rough we should not get MRI in the 
case of application of our numerical method also. 

It is known from the observations that the shapes of 
core collapse supernovae are different. From our simulations 
it follows that MR supernova explosion arises after develop- 
ment of the MRI. The development of the MRI is a stochas- 
tic process and hence the resulting shape of the supernova 
can vary. We may conclude that MR supernova explosion 
mechanism can lead to different shape of the supernova. It 
is important to point out that MR mechanism of supernova 
explosion leads always to asymmetrical outbursts. 

The simulations of the MR supernova explosions for 
the initia l quadrupole-like m agnetic field described in the 
paper by lArdelian et"ai] Jiooai and MR supernova explo- 
sion described in current paper are restricted by the sym- 
metry to the equatorial plane. While in reality this sym- 
metry can be violated due to the MRI, simultaneous pres- 
ence of the initial dipo le and quadrupole -like magnetic 
field llWa uE et al] il992i) ) and in it ial to roidal magnetic field 
iSisnovatvi-Kogan fc Moiseenkd lll992f) ). The violation of 
the symmetry could lead to the kick effect and formation 
of rapidly moving radio pulsars. 

When rotational and magnetic axes do not coincide the 
whole picture of the explosion process is three dimensional. 
Nevertheless, the magnetic field twisting happens always 
around the rotational axis, so we may expect the kick ve- 
locity of the neutron star be strongly correlated with its 
spin direction, also du e to th e anisotropy of the neutrino 
flux iBisnovatvi-Koeanl (^g^)- Simultaneously, because of 
the stochastic nature of MRI, the level of the anisotropy 
should be strongly variable, leading to a large spreading in 



the the neutron star velocities. This prediction of MR ex- 
plosion differs from the models with a powerful neutrino 
convection, where arbitrary direction of the kick velocity is 
expected t Burrow s e^^ (^g^) ) . R ecent analysis of obser- 
vations of pulsars Jjohnston^tal] i2005l) 'l shows, that ro- 
tation and velocity vectors of pulsars are aligned, as it is 
predicted by the MR supernova mechanism. The alignment 
of the vectors can be violated in the case when supernova 
explodes in a binary system. 
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